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1 . SUMMARY 


Durtivg this quarter^ the new "discrete shield" temperature model 
was completed and verified. Modltlcatlona to the J419 low-stress 
configuration were tested experimentally to evaluate effects on growth 
speed. A composite lid and shield configuration combining the low- 
stress features of the J419 with the width-limiting characteristics of 
the J98M3 was fabricated and tested in the N-furnace. Several long 
crystals were grown with width limited to about 3.3 cm and with melt 
replenlslsnent, although the configuration Is not yet optimized for 
steady-state growth. 
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2. INTRODUCTION 


Silicon dendritic web is a single-crystal silicon ribbon 
material which provides substantial advantaf.es for low-cost manufacture 
of solar cells. A signiflcai.t feature of tho orocess is the growth from 
a melt of silicon without constraining dies, resulting in an oriented 
single-crystal ribbon having excellent surface features. In common with 
other more classical processes such as Czochralskl growth. Impurity 
rejection into the melt permits the use of less pure "solar grade" 
starting material without significantly affecting cell performance. A 
unique property of the dendritic web process is the growth of long 
ribbons of controllable width and thickness which not only facilitates 
automation of subsequent processing into solar cells, but also results 
in high material utilization since cutting and polishing are not 
required. 

On the present contract, three broac, areas of work are 
emphasized: 

1 . The development of thermal stress models in order to understand 
the detailed parameters which generate buckling stresses, ''he 
model can then be used to guide the design of Improved low- 
stress web growth configurations for experimental testing. 

2. Experiments to increase our under jtanding of the effects of 
various parameters on the web growth process. 

3. The construction of an experimental web growth machine which 
contains in a single unit all the mechanical and electronic 
features developed previously so that experiments can be carried 
out under tightly controlled conditions. 
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Thus» the principal objective of this work has been to expand our 
knowledge and understanding of both the theoretical and experimental 
aspects of the web growth process to provide a solid base for substan- 
tial Improvements In both area throughput and web crystal quality* 

During this reporting period, a new model for calculating the 
temperature distribution In the growing web was completed. This model 
can handle any number of discrete lids and top shields and thus permits 
much more complex configurations to be accurately represented than was 
possible with the old "lumped shield" model. The model-generated low- 
stress growth configuration was adapted for melt replenishment and width 
control and Installed in the N-furnace for evaluation. 


3. TECHNICAL PROGRESS 


3.1 Improved Heat Transfer Model 

3.1.1 Introduction 

The modeling of the chcnnal stress and buckling of dendritic web 
crystals in fact involves three distinct computer codes: 1) a model 

which calculates the temperature distribution along the growing ribbon; 
2) a model which uses the temperature distribution as input and then 
calculates the thermal stress; and 3) a model which uses the thermal 
stress distribution as input and then calculates the stability of the 
ribbon with respect to buckling. With such a sequential type of 
analysis » it is of obvious importance to have as good as possible a 
model in the first stage. In fact, the second and third stages of 
modeling utilize a general, proprietary finit- element code, WECAN, 
which has been shown to be more than sufficiently accurate for our 
present purposes. The present work concerns the improvement of the 
calculation of the temperature dlstrlbu^:ion in the web which forms the 
basis for the subsequent calculations. 

The model which calculates the temperature distribution along 
the growing silicon ribbon has two distinct routines: a routine which 

calculates the radiative heat exchange between a point on the ribbon and 
its environment and a second routine that integrates the heat conduction 
equation. The problems which occasioned the present work lay not in the 
integration routine, which had been upgraded earlier, but in the 
geometrical factors. 

For several years, the geometry shown in Figure 1 has been used 
to evaluate the radiation Interchange with the ribbon. This model quite 
adequately represented the simpler growth configurations used urtil 
recently, which had only one or two shields above the lid. Currently, 
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Figure 1. Geometrical and thermal parameters of previous model for 
calculating temperature profile In the web. 


ORIGINAL PAGE iS 
OF POOR QUALITY 


growth configurations with multiple, and occasionally widely spaced, 
shields are being used, and representation of such a configuration by an 
Isothermal rectangle is suspect. During the reporting period, we have 
developed a new computer code which allows calculation of the radiative 
Interchange between the web and a number of shields, spacers, etc. so 
that we have greatly Improved the geometric resolution of the model. 

in the previous model, we have neglected other heat transfer 
mechanisms such as gaseous conduction and convection. Furthermore, we 
still consider that the absorbance of the silicon ribbon is equal to the 
emittance, and that the shields, lids, etc. ere black bodies so chat 
multiple reflective transfer is neglected. While such refinements might 
be desired from a theoretical standpoint, the success of the previous 
modeling suggests that they are second order effects. 

3.1.2 Review of Radiative Heat Transfer 

Consider the radiative transfer between a single shield at 
Temperature T^j^ and an element dx of the web at the point at 
temperature T(x) (which is to be determined). Figure 2 illustrates the 
shield subtending an angle from 6 ^^^^ to 6^2 ffom the normal to the 
element dx. The radiative flux density from the shield to the point x 


on the web is 


1 . i2 

y 0 T^ / COS0 d9 


where o is the Stefan-Boltzman constant. The radiative flux density 
from the point x to the shield is « 

I 4 

— eo T(x) j cosP d 6 
^ 0 

il 


where e is the emissivity of silicon web. Since we arc assuming the 
absorbance of the web equals its emittance, the new flux density 
transferring from x is 
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q - Eo [Tfx)^ - T^^](iln 0^2 - gln0j^^)/2 ^ 3 ^ 

The net flux density transfer from any number of elements Is the 
31BB of terms like the one above* The essential fact is that once the 
angulir limits and temperature of each element are determined, q may be 
simply evaluated for substitution in the heat transfer equation^: 

„ dT d -a dT 2q 

^ “ (4) 

where p " density; C ■ specific heat; u • pull velocity; a ■ 318 W/cm.> 
the coefficient of the temperature dependence of thermal conductivity, 
and t web thickness. As before, this equation may be integrated to 
determine the temperature distribution T(x) along the web. Vte need only 
evaluate the geometrl'* factors of equation 3. 

3.1.3 New Geometry 

We examine now the original geometry of Figure 1 to see how it 
can be systematized for extending Its number of elements. A possible 
systematic arrangement is Illustrated in Figure 3; by numbering the 
ambient regions in sequence with the lid and shield i, they may be 
treated in the same fashion as the other elements. The melt region and 
the crucible wall are assumed to be at the same temperature and are 
assigned the number "O'*. In this case, the ambient region above the top 
shield would be numbered element 5. The lid is divided into two 
elements because of the bevel. Since the heat transfer depends only on 
angular subtense (reflections are neglected), the beveled lid Is 
equivalent to a stepped lid as Indicated by the dotted lines. In the 
previous gecxnetry, the whole lid was isothermal, but now, if desired, 
different temperatures can be assigned to the top and bottom elements of 
the lid. Ii any case, replacing the bevel with a step makes all the 
elements rectangular in shape. Since the geometries of interest have 
the back edges of the elements lining up, the ambient spaces can be 
considered to be comprised of rectangles such as elem«:.nt 3. Thus, the 
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Figure 3. Shading limits for viewing the furnace elements from the web. 
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geometry is completely determined by the Input of the upper corner of 
each element Indicated by the large, black points* 

The four numbered points on the web are the boundaries of 
regions with different views; for example, for points between 3 and 4, 
all elements are In view, while above 4, the ambient region 3 Is 
Invisible. These points along the x-axls can be systematically 
determined in the following manner. First, the lower "shade" of each 
element Is fouiid by drawing lines through the inner upper corner of each 
element below It. The highest intersection of these lines on the x-axis 
is the lower shade. For example, the point 1 is the lower shade of the 
4th element; element 4 is Invisible for points on the web below point 1. 

After all the lower shades are found (some elements may not have 
any) , then the upper shades are found. The lower inner corner of an 
element Is connected to the lower Inner c<'rner of each element above 
It. The lowest Intersection point of each of these lines with the 
x-axls Is the upper shade of the element in question. For example, 
element 2 does not have an upper shade, while element 3's upper shade Is 
the point 4 on the web; element 3 Is Invisible above point 4 on the web. 

After all the upper and lower shades are found, these points are 
ordered on the x-axls. The region between two such points is assigned a 
vector, the l^h component of which Is either one or zero depending on 
whether or not the 1th element Is visible to points In this region. 

Figure 4 illustrates a possible geometric configuration for the 
new model. Numerical integration of the heat conduction equation 
proceeds up the silicon web from the melt surface. By knowing in which 
region between upper and lower shades the point of Integration lies, we 
know which elements are visible to it. For the point illustrated, 
elements (the melt), 4, 6, 7, 8, aid 15 (ambient) are visible. There 
Is one additional piece of information needed before the Integration can 
proceed; the e lent opposite the point of Integration, element 7 In 
this case. For the visible elements numbered 7 or less, the lower Inner 
corners are used to limit the viewing angles while tor elements 7 or 


10 





ORIGINAL PAGE 13 
OF POOR QUALITY 


greater, the upper inner corners are used as Illustrated. The viewing 
regions are labeled 1 through 6; these determine the proper coefficient 
q of the heat conduction equation during integration. 

Once the geometrical heat transfer coefficients have been calcu- 
lated, equation 4 can be numerically Integrated to yield the temperature 
distribution along the length of the web crystal. A more detailed 
description of the computer code will be presented in a future report. 

3.1.4 Verification of the New Model 

One of the modeling runs of the J419 configuration. Case 9-llC, 
was used to verify the new radiation transfer model. The first test was 
to run the identical configuration, i.e., a lumped shield model as in 
Figure 1, to verify the operation of the routine which evaluates the 
shade points, etc. The results were excellent; the new program gave the 
identical temperature profile as the older program. A second geometry 
was then evaluated in which both the lid and the shield "block" were 
bei?eled. This is a relatively small change in the geometry, but one 
that could not be accommodated by the older model. It would be expected 
that the results be similar, but not identical, and that was Indeed the 
case. 

Figures 5 and 6 illustrate the two cases.* However, Instead of 
presenting the temperature profiles themselves, the figures show the 
second derivative of temperature times thermal expansion coefficient, 
i.e., (aT)". This parameter was chosen because It accentuates the 
cur\ature of the temperature which can basically be considered as the 
cause of thermal stress. Examination of the two figures shows that 
indeed the two profiles are almost the same, although there are some 
slight differences due to the more complicated profile. 


♦These figures were generated by a computer plotting routine which has 
been added to the temperature distribution model. Not only is (aT)" 
plotted, but also the lid and shield geometry is represented (heavy 
lines in the figures). 


12 


-AXIS 


ORIGINAL PAGE 13 
OF POOR QUALITY 



Figure 5. Model geometry and (aT)" for run ^-IIC. 
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Figure 6. Model geometry and (aT)" for run 9-llCA (beveled elements). 
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3 . 2 Growth Experiments 

During the early part of this reporting period, growth runs In 
the new experimental facility (N~Furnace) were dedicated to achieving 
the proper furnace adjustments for continuous feeding of a standard 
width-limiting J98M3A lid and shield assembly combined with an elongated 
crucible. Previous runs had Indicated that when the adjustable end 
shield position was proper for melting the feed pellets, then the melt 
profile was too dipped (colder in the center of the growth slot than at 
the ends), which led to growth instability, pull outs, and extra 
dendrite generation. Increasing the diameter of the outboard end holes 
In the lid flattened out the melt ptoflle and greatly improved the 
general growth. 

Although the melt level control system components all 
functioned. It became apparent that the system as a whole was not 
properly tuned to the growth requirements. Several changes were made In 
the control loop which brought the response time and loop gain Into 
better agreement with the control requirements, and excellent results 
were then obtained with actual web growth. In one run, several crystals 
were grown under automatic level control. During the whole course of 
the run, all the crystals maintained a constant thickness within about 
4 pm (at one constant growth speed) and even the hold temperature did 
not change by more than a few microvolts. The extreme stability of the 
growth conditions would Indicate that the level control system was 
operating even better than required! 

The modifications made to the N-furnace melt level control 
system to better meet growth requirements were three In number. First, 
the time constant of the control loop was made adjustable so that 
Integration periods of 1, 20, 100, and 400 seconds could be used. The 
faster response times are used for the initial set up of the equipment, 
but the longer times are needed to avoid unnecessary fluctuations In the 
final feed rates. With the longer Integration tiroes, a higher loop gain 
can also be used and the response gain of the sensing circuit was 
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increased by a factor of about ten. As a result, the control 
capabilities of the system are easily of the order of 100 pm or better 
in melt height sensing. Finally, a selectable limiting speed was placed 
on the pellet feed mocor in order to prevent overfeeding during 
transient conditions. 

Several modifications of the computer model-generated J419 
configuration were tested in the J-furnace in an effort to enhance the 
growth speed of this wide-growth, low-stress design. In one series of 
experiments, the J419 was modified to a more open geometry by beveling 
the top lid and using a different shield stack order. The expected speed 
enhancement occurred, as did a change in the residual stress character. 
Whereas the standard J419 configuration produced crystals with negative 
residual stress, the modified version gave material with positive 
residual stress. In both cases, the magnitude could be either very small 
or moderate, depending on other growth parameters. The higher speed 
modifications also would grow moderately wide material, up to 40 mm, but 
tended to deform at slightly narrower widths than the unmodified J419. 

The earlier modeling results suggested that the residual stress 
in some of the higher speed versions of the configuration could be 
reduced by reducing the heat loss from the web in the Interface region. 
One feature of a lid design which would accomplish this would be a bevel 
on the bottom edge of the slot in the lid. Thermally, this would allow 
the web to "see" the hot crucible cavity for a greater portion of its 
length so that the heat loss near the growth front would be decreased. 

Four runs (J439-J442) were made in which the top lid was either 
unbeveled or had a 3 mm bevel; the bottom lid had either a 1.3 mm or a 3 
mm bevel. It was found that the bevel on the bottom of the lid had more 
effect on speed than on residual stress, while the bevel on the top lid 
affected Che stress more chan the speed. Although the previous thermal 
model does not have enough resolution to evaluate these effects, the new 
model should easily accommodate these geometry changes and will be used 
to provide the required insight into the behavior of the system. 


While modifications to the J4I9 configuration were being tested 
in the J-furnace» experiments in the WA-furnace were directed to testing 
the compatibility of the J419 configuration with the elongated crucible 
which would be required for growth at constant melt level* Adjustable 
end shields were installed. Melt probe data showed that the lateral 
coll position and the end shield height affect the melt profile 
differently and that adjustment of both parameters is necessary to 
obtain a flat, symmetric temperature distribution In the melt. However, 
when the proper melt temperature profile was achieved, it was 
established that the low-stress characteristics of the J419 
configuration were maintained with the elongated crucible configuration. 

At this point it seemed appropriate to combine the wlde-growth, 
low-stress capabilities of the J419 configuration with the width-limiting 
behavior of the J98M3A configuration for demonstrating extended growth 
runs. Lids and shields were fabricated which combined these features in 
a hybrid configuration now known as the J435. Four runs were made with 
this configuration in the J-Furnace to evaluate the width-limiting 
capabilities of the design, and the results confirmed a steady growth 
width between 28 ram and 33 mm, depending on the furnace parameters. 

Feed holes and laser holes were then added to the J435 lids and 
shields and the hardware Installed in the N-Furnace with an elongated 
crucible for melt replenishment. Preliminary experiments focused on 
optimizing the end shield adjustment and the work coll position so that 
the melt profile is comoatlble with both steady-state width and melt 
replenishment . 

Several growth runs were made both with and without melt 
replenishment. In one run, a 4.8 meter crystal was grown with the width 
held at 3.1 to 3.3 cm for 3.4 meters. This crystal was grown without 
replenishment due to clogging of the feed hole with oxide. Several 
shorter crystals were grown with replenishment, trith thickness constant 
to within a few microns and width limited in the 3.1 to 3*3 cm range. 
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We found tliat clogging of the feed hole with oxide can be 
prevented by constantly feeding pellets at a slow rite throughout the 
run, whether a crystal Is growing or not. We will also test some 
modifications In the feed hole geometry to alleviate this problem. 


4. CONCLUSIONS 


The new temperature model will allow complex lid and shield 
configurations to be evaluated for stress-inducing characteristics with 
the computer* The greatly Improved resolution also permits minor design 
variations to be tested so that only the most promising designs need be 
fabricated into hardware for experimental evaluation. Thus, the 
experimental effort can be concentrated on designs with high probability 
of success and on refinements which can only be tested empirically. 

The J435 configuration, which combines the stress-reducing 
characteristics of the J419 with the width-limiting capabilities of the 
J98M3, has been tested in the N-furnace and appears capable of producing 
long crystals about 3.3 cm in width. Further refinements are needed to 
better control oxide deposition, and the optimum melt level for steady- 
state growth is yet to be determined. 
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5. PLANS AND FUTURE WORK 


The new temperature model will be used to generate new lower 
streso growth configurations. Experimental work will continue to 
optimise the J43S configuration for steady~state growth. 
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6. NEW TECHNOLOGY 


No new technology Is reportable for the period covered. 
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